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Abstract 

A first flavour-tagged measurement of the time-dependent CP- violating asymmetry 
in Bg — > (f)(f) decays is presented. In this decay channel, the CP-violating weak phase 
arises due to CP violation in the interference between B®-B® mixing and the b — > sss 
gluonic penguin decay amplitude. Using a sample of pp collision data corresponding 
to an integrated luminosity of 1.0 fb _1 and collected at a centre-of-mass energy of 
7 TeV with the LHCb detector, 880 B® — > <f)(j) signal decays are obtained. The 
CP-violating phase is measured to be in the interval [—2.46, —0.76] rad at 68% 
confidence level. The p-value of the Standard Model prediction is 16%. 
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The — > (fxf) decay is forbidden at tree level in the Standard Model (SM) and proceeds 
via a gluonic b — > sss penguin process. Hence, this channel provides an excellent probe of 
new heavy particles entering the penguin quantum loops [1-3]. Generally, CP violation in 
the SM is governed by a single phase in the Cabibbo-Kobayashi-Maskawa quark mixing 
matrix [4] . The interference between the Bg-B° s oscillation and decay amplitudes leads to a 
CP asymmetry in the decay time distributions of B® and B® mesons, which is characterised 
by a CP-violating weak phase. The SM predicts this phase to be small. Due to different 
decay amplitudes the actual value is dependent on the B° s decay channel. For £?° — > J/ifxf), 
which proceeds via a b — > ccs transition, the SM prediction of the weak phase is given 
by -2 arg (-V ts V* b /V cs V* b ) = -0.0364 ± 0.0016 rad [5]. The LHCb collaboration recently 
measured the weak phase in this decay to be 0.068 ± 0.091(stat) ± O.Oll(syst) rad [6], 
which is consistent with the SM and places stringent constraints on CP violation in B®-B° s 
oscillations [7]. In the SM, the phase in the 5° — >■ 00 decay, S , is expected to be close 
to zero due to a cancellation of the phases arising from B®-B° s oscillations and decay [8]. 
Calculations using QCD factorization provide an upper limit of 0.02 rad for |0 S | [1-3]. 

In this Letter, we present the first measurement of the CP- violating phase in B® — > (fxf) 
decays. Charge conjugate states are implied. The result is based on pp collision data cor- 
responding to an integrated luminosity of 1.0 fb" 1 and collected by the LHCb experiment 
in 2011 at a centre-of-mass energy of 7 TeV. This data sample was previously used for 
a time-integrated measurement of the polarisation amplitudes and triple product asym- 
metries in the same decay mode [9]. The analysis reported here improves the selection 
efficiency, measures the B° s decay time and identifies the flavour of the B° s meson at pro- 
duction. This allows a study of CP violation in the interference between mixing and decay 
to be performed. 

The LHCb detector is a forward spectrometer at the Large Hadron Collider covering 
the pseudorapidity range 2 < r\ < 5 and is described in detail in Ref. [10]. Events 
are selected by a hardware trigger, which selects hadron or muon candidates with high 
transverse energy or momentum (pt), followed by a two stage software trigger [11]. In the 
software trigger, B° s — > 00 candidates are selected either by identifying events containing a 
pair of oppositely charged kaons with an invariant mass close to that of the meson or by 
using a topological 6-hadron trigger. In the simulation, pp collisions are generated using 
Pythia 6.4 [12], with a specific LHCb configuration [13]. Decays of hadronic particles 
are described by EvtGen [14], the detector response is implemented using the Geant4 
toolkit [15] as described in Ref. [16]. 

The B® — > <f)<f> decays are reconstructed by combining two meson candidates that 
decay into the K + K~ final state. Kaon candidates are required to have p-r > 0.5 GeV/c, 
impact parameter (IP) x 2 larger than 16 with respect to the primary vertex (PV), where 
the IP x 2 is defined as the difference between the x 2 °f the PV reconstructed with and 
without the considered track. Candidates must also be identified as kaons using the ring- 
imaging Cherenkov detectors [17], by requiring that the difference in the global likelihood 
between the kaon and pion mass hypotheses (AIylCktt = ln£^- — lnC w ) is larger than 
—5. Both meson candidates must have a reconstructed mass, ttikk, of the kaon pair 
within 20 MeV/c 2 of the known mass of the meson, a transverse momentum (p^) larger 
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than 0.9 GeV/c and a product p^ 1 ■ pi^ > 2 GeV 2 /c 2 . The \ 2 P er degree of freedom (ndf) 
of the vertex fit for both meson candidates and the B® candidate is required to be 
smaller than 25. Using the above criteria, 17575 candidates are selected in the invariant 
four-kaon mass range 5100 < ttikkkk < 5600 MeV/c 2 . 

A boosted decision tree (BDT) [18] is used to separate signal from background. The 
six observables used as input to the BDT are: Pt> V an d x 2 /ndf of the vertex fit for the B® 
candidate and the cosine of the angle between the B° s momentum and the direction of flight 
from the closest primary vertex to the decay vertex, in addition to the smallest pr and 
the largest track x 2 /ndf of the kaon tracks. The BDT is trained using simulated B° s — > 00 
signal events and background from the data where at least one of the candidates has 
invariant mass in the range 20 < \mxK — ^<j>\ < 25 MeV/c 2 . 

The sPlot technique [19,20] is used to assign a signal weight to each B® — > 00 candidate 
in order to separate signal from background, employing the four-kaon mass as the discrim- 
inating variable. The sensitivity to S is optimised taking in to account the signal purity 
and the flavour tagging performance. The final selection of B° s — > 00 candidates based 
on this optimisation is required to have a BDT output larger than 0.1, A\nC Kn > —3 for 
each kaon and \rriKK — m $\ < 15 MeV/c 2 for each candidate. 

In total, 1182 B® — > 00 candidates are selected. Figure [I] shows the four-kaon in- 
variant mass distribution for the selected events. Using an unbinned extended maximum 
likelihood fit, a signal yield of 880 ± 31 events is obtained. In this fit, the B° s — > 00 signal 
component is modelled by two Gaussian functions with a common mean. The relative 
fraction and width of the second Gaussian component are fixed to values obtained from 
simulation. The width of the first Gaussian component is 12.9 ± 0.5 MeV/c 2 , in agree- 
ment with the expectation from simulation. Combinatorial background is modelled using 
an exponential function which is allowed to vary in the fit. Contributions from specific 
backgrounds such as B° — > <pK*°, where K*° — > K + tt~ , are found to be negligible. 

An unbinned maximum likelihood fit is performed to the decay time, t, and the three 
helicity angles, Q = {9i, 9 2 , $} [9], of the selected B° s — > 00 candidates, each of which is re- 
assigned a signal sPlot weight based on the four-kaon invariant mass, ttikkkk [19,20]. The 
probability density function (PDF) consists of signal components, which include detector 
resolution and acceptance effects, and are factorised into separate terms for the decay 
time and the angular observables. 

The B° s decay into the K + K~K + K~ final state can proceed via combinations of 
intermediate vector (0) and scalar (/ (980)) resonances and scalar non-resonant K + K~ 
pairs. Thus the total decay amplitude is a coherent sum of P-wave (vector- vector) , S'-wave 
(vector-scalar) and S'S'-wave (scalar-scalar) contributions. The differential decay rate of 
the decay time and helicity angles is described by a sum of 15 terms, corresponding to 
five polarisation amplitudes and their interference terms, 



The angular functions fi(fl) for the P-wave terms are derived in Ref. [21] and the helicity 
angles of the two mesons are randomly assigned to 9\ and 82- The time-dependent 
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Figure 1: Invariant K + K~ K + K~ mass distribution for selected B® — > <f><j) candidates. The 
total fit (solid line) consists of a double Gaussian signal component together with an exponential 
background (dotted line). 

functions Ki(t) can be written as 

Ki(t) = Nie~ Tst [ci cos(Am s t) + d { sin(Am s t) + a* cosh(§Ar s t) + b { sinh(f Aiy)], (2) 

where Ar s = Tl — Th is the decay width difference between the light (L) and heavy 
(H) B® mass eigenstates, T s is the average decay width, T s = (i\ + Th)/2, and Am s is 
the B®-B° s oscillation frequency. The coefficients Ni, a>i, bi, Ci and di can be expressed in 
terms of (p s and the magnitudes, \Ai\, and phases, 8i, of the five polarisation amplitudes 
at t — 0. The three P-wave amplitudes, denoted by Ao,Aii,Aj_, are normalised such 
that | 1 2 + |^4||| 2 + |^4-i_| 2 — 1, with the strong phases 5± and 62 defined as 5± = 5± — 5\\ 
and S2 = S± — 5o. The S and SS-w&ve amplitudes and their corresponding phases are 
denoted by As, Ass and 5s, Sss, respectively. For a B® meson produced at t = 0, the 
coefficients in Eq.[2]and the angular functions fi(Q\, 9 2 , $) are given in TableJTJ where 5 2) i = 
82 — 8\. Assuming that CP violation in mixing and direct CP violation are negligible, the 
differential distribution for a B° s meson is obtained by changing the sign of the coefficients 
Ci and di. The PDF is invariant under the transformation (</>„, Ar s , 5\\, 5±, 8s, Sss) ~^ 
(n — (f) s , — AT S , —8n, ir — 8±, —8s, —8ss)- This two-fold ambiguity is resolved in the fit as the 
B° s average decay width and decay width difference are constrained to the values measured 
in B° s J/ip(f) decays, T s = 0.663 ± 0.008 ps" 1 and Ar s = 0.100 ± 0.017 ps" 1 with a 
correlation coefficient p(Ar s , T s ) = —0.39 [6]. Similarly, the B® oscillation frequency Am s 
is constrained to the value Am s = 17.73 ± 0.05 ps -1 [22]. 

A correction factor is multiplied to the interference terms in Table [l] between the 
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Table 1: Coefficients of the time-dependent terms and angular functions defined in Eqs. [I] and 
[2] Amplitudes are defined at t = 0. 
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P and S'-wave (and the P and SS-w&ve) contributions to account for the finite itikk 
mass window considered in the amplitude integration. This factor is calculated from the 
interference between the different tukk lineshapes of the vector and scalar contributions. 

The acceptance as a function of the helicity angles is not completely uniform due to 
the forward geometry of the detector and the momentum cuts placed to the final state 
particles. A three-dimensional acceptance function is determined using simulation. The 
acceptance factors are included in the fit as a normalisation of the PDF for each of the 
angular terms. The acceptance function varies by less than 20% across the phase-space. 

The event reconstruction, trigger and offline selections introduce a decay time depen- 
dent acceptance. In particular for short decay times, the acceptance vanishes due to the 
trigger, which requires tracks with significant displacement from any PV. The decay time 
acceptance is determined using simulation and incorporated by multiplying the signal 
PDF with a binned acceptance histogram. 

The parameters of a double Gaussian function used to model the decay time resolution 
are determined from simulation studies. A single Gaussian function with a resolution of 
40 fs is found to have a similar effect on physics parameters and is applied to the data fit. 

The (f) s measurement requires that the meson flavour is tagged as either a B° s or B° s 
meson at production. To achieve this, both the opposite side (OS) and same side kaon 
(SSK) flavour tagging methods are used [23,24]. In OS tagging the 6-quark hadron 
produced in association with the signal 6-quark is exploited through the charge of a muon 
or electron produced in semileptonic decays, the charge of a kaon from a subsequent 
charmed hadron decay, and the momentum- weighted charge of all tracks in an inclusively 
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reconstructed decay vertex. The SSK tagging makes use of kaons formed from the s-quark 
produced in association with the B° s meson. The kaon charge identifies the flavour of the 
signal Bg meson. 

The event-by-event mistag is the probability that the decision of the tagging algorithm 
is incorrect and is determined by a neural network trained on simulated events and cali- 
brated with control samples [23] . The value of the event-by-event mistag is used in the fit 
as an observable and the uncertainties on the calibration parameters are propagated to 
the statistical uncertainties of the physics parameters, following the procedure described 
in Ref. [6]. For events tagged by both the OS and SSK methods, a combined tagging 
decision is made. The total tagging power is EtagD 2 = (3.29 ± 0.48)%, with a tagging 
efficiency of e tag = (49.7 ± 5.0)% and a dilution V = (1 — 2oj) where u is the average 
mistag probability. Untagged events are included in the analysis as they increase the 
sensitivity to (f) s through the bi terms in Eq. [2] 

The total S'-wave fraction is determined to be (1.6 where the double S-wave 

contribution Ass is set to zero, since the fit shows little sensitivity to Ass- As a cross-check 
a sideband subtracted fit to the two-dimensional mass, m^K, for both kaon pairs, with 
background subtracted using sidebands, yields a consistent S'-wave fraction of (2.1±1.2)%. 

The results of the fit for the main observables are shown in Table [2] Figure [2] shows 
the distributions for the decay time and helicity angles with the projections for the best 
fit PDF overlaid. The likelihood profile for the CP- violating weak phase <p s , shown in 
Fig. |3j is not parabolic. To obtain a confidence level a correction is applied due to a 
small under- coverage of the likelihood profile using the method described in Ref. [25]. 
Including systematic uncertainties (discussed below) and assuming the values of nuisance 
parameters observed in data, an interval of [—2.46, —0.76] rad at 68% confidence level is 
obtained for S . The polarisation amplitudes and phases, shown in Table |2j differ from 
those reported in Ref. [9] as S is not constrained to zero. 

The uncertainties related to the calibration of the tagging and the assumed values of 
T s , Ar s and Am s are absorbed into the statistical uncertainty, described above. System- 
atic uncertainties are determined and the sum in quadrature of all sources is reported in 
Table [2] for each observable. To check that the background is properly accounted for, an 
additional fit is performed where the angular and time distributions are parameterised us- 

Table 2: Fit results with statistical and systematic uncertainties. A 68% statistical confidence 
interval is quoted for <f> s . Amplitudes are defined at t = 0. 



Parameter 



Value 



fstat. 



^syst. 



fl [rad] (68 % CL) 



[-2.37,-0.92] 
0.033 
0.046 



0.22 
0.017 
0.018 
0.009 
0.12 
0.10 
0.33 




0.329 
0.358 
0.016 

2.19 
-1.47 

0.65 



+0.024 
-0.012 



0~i [rad] 
62 [rad] 
5s [rad] 



0.44 
0.48 

+0.89 
-1.65 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 

cos 6, cos 9 2 

Figure 2: One-dimensional projections of the B® — > (fxf) fit for (a) decay time, (b) helicity 
angle $ and the cosine of the helicity angles (c) 0\ and (d) 82- The data are marked as points, 
while the solid lines represent the projections of the best fit. The CP-even P-wave, the CP-odd 
P-wave and 5-wave components are shown by the long dashed, short dashed and dotted lines, 
respectively. 



ing the B® mass sidebands. This gives results in agreement with those presented here and 
no further systematic uncertainty is assigned. The uncertainty due to the modelling of 
the S"-wave component is evaluated by allowing the SS- wave component to vary in the fit. 
The difference between the two fits leads to the dominant uncertainty on <p s of 0.20 rad. 
The decay time acceptance is explicitly parameterised and, by taking the difference in 
the observed value of <f> a , a systematic uncertainty of 0.09 rad is determined. Possible 
differences in the simulated decay time resolution compared to the data are studied by 
varying the resolution according to the discrepancies observed in the B® — > J/ip<fr anal- 
ysis [6]. This leads to a systematic uncertainty of 0.01 rad for <p s . The distributions of 
maximum and xV n df of the final state tracks and the and V of the B° s candidate are 
reweighted to better match the data. From this, the angular acceptance is recalculated, 
leading to small changes in the results (0.02 rad for <p s ), which are assigned as systematic 
uncertainty. Biases in the fit method are studied using simulated pseudo-experiments 
that lead to an uncertainty of 0.02 rad for S . Further small systematic uncertainties 
(0.02 rad for <f> s ) are due to the limited number of events in the simulation sample used 
for the determination of the angular acceptance and to the choice of a single versus a dou- 
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-2 2 

$ t [rad] 

Figure 3: Negative Aln likelihood scan of (fi s . Only the statistical uncertainty is included. 



ble Gaussian function for the mass PDF, which is used to assign the signal weights. The 
total systematic uncertainty on <fr s is 0.22 rad, significantly smaller than the statistical 
uncertainty. 

In summary, we present the first study of CP violation in the decay time distribution 
of hadronic B® — > cfxp decays. The CP-violating phase, <p s , is restricted to the interval 
of [—2.46, —0.76] rad at 68% C.L. The p-value of the Standard Model prediction [8] is 
16%, taking the values of the strong phases and polarisation amplitudes observed in data. 
The precision of the (p s measurement is dominated by the statistical uncertainty and is 
expected to improve with larger LHCb data sets. 
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